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C/3 1 Abstract. 

>v 

Motivated by many observations of anomalies in condensed matter systems, we 
consider a new fundamental Hamiltonian in which condensed matter and nuclear 
systems are described initially on the same footing. Since it may be possible that 
the lattice will respond to the mass change associated with a excited nuclear state, we 
adopt a relativistic description throughout based on a many-particle Dirac formalism. 
This approach has not been used in the past, perhaps due to the difficulty in separating 
the center of mass and relative degrees of freedom of the nuclear system, or perhaps 
due to an absence of applications for such a model. In response to some recent ideas 
about how to think about the center of mass and relative separation, we obtained 
from the Dirac model a new fundamental Hamiltonian in which the lattice couples to 
different states within the composite nuclei within the lattice. In this description the 
different nuclear states have different mass energies and kinetic energies, as we had 
expected. In addition there appear new terms which provide for nuclear excitation as 
a result of coupling to the composite momentum. This new effect comes about because 
of changes in the composite nuclear state as a result of the dynamical Lorentz boost 
in the lattice. 



1. Introduction 

From a condensed matter viewpoint, a solid is made of nuclei and electrons, where the 
nuclei can for the most part be treated as point particles pQ. In cases where a more 
sophisticated description is needed, the basic description is augmented with nuclear 
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spin, magnetic moments, and electric quadrupole moments [2j. Essentially the only 
place that excited states show up is in Mossbauer studies where they are required 
in order to describe the absorption or emission of a gamma [3]. The picture that 
results is both wonderfully simple and very rich; simple, in that the Born-Oppenheimer 
separation allows us to reduce the problem to electronic bands and phonon modes; and 
rich, as the models that result describe a wide range of basic, subtle, and occasionally 
unexpected physical phenomena. This basic approach to condensed matter physics has 
been sufficiently successful over the years that it would require some rather dramatic 
new experimental result before we might be motivated to revise it in any significant 
way. 

In recent years there have been claims of experimental results which at a 
fundamental level seem not to be consistent with this basic viewpoint of condensed 
matter physics. In the Fleischmann-Pons excess heat effect [U, [5], a great deal of 
energy is generated which is thought to have a nuclear origin (since there are no 
commensurate chemical byproducts, and since 4 He is observed as a possible product in 
amounts proportional to the energy produced [6j[7J[8]), without commensurate energetic 
nuclear products [9]. In these experiments it almost seems as if the solid is taking 
up an MeV quantum; if so, then this constitutes an effect which seems very hard to 
understand within our current condensed matter framework. Given that such an effect 
seems impossible to contemplate within modern condensed matter physics (and also 
within modern nuclear physics), a natural reaction has been to go with the existing 
picture (supported by a very large body of experimental results and a consistent and 
mature theory), and to reject the Fleischmann-Pons experiment as simply being in error 

During the past two decades and more this de facto solution has been adopted 
generally, and it has worked surprisingly well. Science has advanced substantially; 
there are now even more experimental and theoretical results which support the modern 
condensed matter viewpoint; and those who pursue anomalies such as the Fleischmann- 
Pons experiment are isolated from science and ignored [11] . 

Meanwhile, another experimental result has been put forth which challenges our 
modern view of condensed matter physics. Karabut has studied a variety of anomalies in 
high-current density glow discharge experiments, and in the course of the work noticed 
that collimated X-ray emission occurred in powerful bursts normal to metal coin-shaped 
samples that served as cathodes |12j [131 El US]- Although anomalous emission effects 
are seen when the discharge is on, the powerful bursts of collimated X-rays are observed 
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on the order of a millisecond after the discharge is turned off. A current spike occurs 
when the discharge is shut off, which we might imagine causes vibrations in the sample. 
It seems as if the vibrational energy is being communicated somehow to produce in- 
phase electronic or nuclear excitation at X-ray energies, which then produces collimated 
X-ray emission through a phased array emission effect. We note that a related effect 
involving the collimated emission of gamma-rays in beamlets was reported earlier by 
Gozzi [16]. Needless to say, such an effect has no place in modern condensed matter 
physics. 

In experiments performed by the Piantelli group, hydrogen is absorbed in nickel 
samples at elevated temperature, resulting in a thermal effect (consistent with energy 
generation) [17], [H], low- level nuclear effects (gamma and neutron emission |19].|20j). 
and the appearance of new elements [21] . This latter effect (appearance of new elements) 
in these experiments is not a low- level effect. Once again, such effects are not predicted 
in modern condensed matter physics. 

These experimental results, and many others, have motivated us to explore new 
models that might be relevant. A major issue that we have been interested in is the 
possibility of coherent energy exchange between quantum systems with mismatched 
characteristic energies, which we considered to be the biggest theoretical problem 
associated with the anomalies. Coherent energy exchange between mismatched quantum 
systems occurs in high harmonic generation [22], so we know that it is possible in 
principle. However, there seems to be no analog to Corkum's mechanism [23j,[24j 
present in the condensed matter system. A lesser version of the effect is known within 
the multiphoton regime of the spin-boson model, which is used to model basic linear 
interactions of two- level systems with an oscillator [25], [26], [27]. We found that if the 
two-level system is augmented with loss, the coherent energy exchange rate is increased 
dramatically. This is due to the fact that destructive interference limits the rate at 
which coherent energy exchange occurs in the spin-boson model, so augmenting the 
model with a mechanism that removes this destructive interference would be expected 
to improve coherent energy exchange rates |28j.[29].|30].|31j. 

Coherent energy exchange in these models works best when the coupling between 
the two-level transition (representing electronic and nuclear transitions) and oscillator 
(representing a vibrational mode) is strong. We studied a further generalization of the 
lossy spin-boson model in which two transitions are coupled to an oscillator, where one 
is strongly coupled and one is weakly coupled [32]. We found that the strongly coupled 
system could assist coherent energy exchange for the weakly coupled system. The 
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model that resulted appeared to us to be very closely related to excess heat production 
in the Fleischmann-Pons experiment, assuming that the mechanism involved D2/ 4 He 
transitions that were weakly coupled to a phonon mode (weakly coupled due to the 
Coulomb repulsion between the deuterons), and that a strongly coupled transition were 
also present. The big problem in this kind of model ends up being the identification 
of the strongly coupled transition. Finding an appropriate strongly coupled transition 
with sufficiently strong coupling to do the job seems problematic within the approach 

After analyzing many candidate transitions, we came to the conclusion that 
there were no physical transitions which could serve as the strongly coupled two-level 
transition within the model. We were optimistic in our writing about the possibility 
that systems described by three-level systems (or iV-level systems) would be able to do 
the job. After putting in a great deal of work on analyzing the strongly coupled three- 
level system, it seemed once again that the coupling was simply not strong enough to 
make a connection with experiment. This conclusion was supported by spectral data 
from Karabut [33], which seemed to be qualitatively consistent with the approach and 
models, but which would require much stronger coupling to explain. 

All of this has led us to the conclusion that we are going to need a new kind of model 
in order to account for the experimental results. To obtain coupling sufficiently strong to 
be consistent with the Karabut experiment, we require very strong interactions that are 
on the general order of what would occur in a nuclear configuration mixing calculation. 
Yet there is no reason to expect that nuclear configuration interactions can couple to 
a phonon mode. In our earlier efforts to describe such an effect, we concluded that 
the internal nuclear degrees of freedom associated with configuration mixing separated 
cleanly from the vibrational degrees of freedom. For many years there has seemed to be 
no viable solution within the general approach, which has been very discouraging. 

The intuitive picture that has emerged over the past few years of thinking about 
the problem is that the different excited states of the nucleus have different masses, 
and under appropriate conditions it may be possible for the nucleus to notice the mass 
differences of the different configurations. This could provide the physical basis for 
phonon exchange in association with configuration mixing. To describe such an effect, 
we need to develop a description of the associated coupling, which seems not to be 
available in the literature. One approach is to begin with a relativistic model for the 
nucleons, and then reduce it in some way to obtain a low momentum approximation in 
which the associated mass effects are retained. 
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The issues under discussion fit within the generic heading of relativistic quantum 
mechanics, which before 1950 would likely have implied the two-body Dirac equation 
as a starting place. However, the need for a manifestly covariant relativistic quantum 
theory more generally led to the development of modern quantum field theory, which 
could in principle be used for problems of interest to us. Field theory is much more 
complicated than relativistic quantum mechanics, so we would prefer a simpler model 
derived from relativistic quantum mechanics if possible. In this day and age, there are 
many relativistic quantum models that have been derived from field theory (such as the 
Bethe-Salpeter equation, as well as others [35J). 

The separation of relative and center of mass degrees of freedom that is possible 
in the nonrelativistic problem does not extend to the relativistic problem, which 
complicates things. In a recent manuscript [36] we obtained a result with the two- 
body Dirac model that seemed to suggest it was possible to arrange for a separation of 
the relative and center of mass contributions to the energy in a simple way. We originally 
made use of this approach in the this paper to obtain an approximate model for the 
center of mass and relative dynamics discussed in this paper. However, in responding 
to the reviewers' comments on the paper, it became clear that there were weaknesses in 
the approach outline, and our interpretation of the results have changed. We were able 
to revise this paper in the galley proof stage in order to take into account the improved 
point of view. The model that results (which is essentially the same in both cases) can 
then be used directly to develop a new Hamiltonian for nuclei in a lattice that includes 
the coupling consistent with a many-particle Dirac formulation. 

Interestingly, the model that results seems to include a relativistic effect which 
provides a direct coupling between the lattice motion and excitations in the nucleus. 
The resulting model appears to be much more closely connected with our earlier models 
than we had expected, which provides the motivation to explore the model further in 
the future. 

2. Center of mass and relative contributions to the energy 

We begin with a description of the free nucleus in terms of Dirac particles within the 
context of a many-particle Dirac Hamiltonian 



Protons and neutrons are composite particles with internal quark structure, and one 
might criticize the use of a Dirac point-particle description for composites in this case; 




(1) 
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however, for our purposes it seems the simplest place to start. The Dirac a and j3 
matrices are 

The interaction between two nucleons appears here as Vjk{vk — r?); we assume that this 
includes strong force and electromagnetic interactions. We assume the $ is an exact 
solution to the time-independent equation 



E<$> 



E a i • °Pj + Pi M 3 c2 + E V 3k( r k - Tj) 

j 3<k 



$. (3) 



2.1. Center of mass and relative coordinates 
The classical center of mass coordinate satisfy 

MR = 5>,t,. (4) 

j 

with 

M = (5) 
j 

The relative position coordinates are 

= r, - R. (6) 
The total classical momentum is 

P = Ep; (7) 

j 

and the relative momenta are 

*i = Pi ~ ^P- (8) 

One of the relative position operators is redundant, since the sum of all relative position 
operators is zero; similarly one of the relative momentum operators is redundant. 
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2.2. Mass parameter 

The eigenvalue can be expressed in terms of relative and center of mass matrix elements 
according to 



E 



E 



m, 



M 3 



cP 



+ ($ £ <*i • c ^ + E /W + £ J^tafc - * ), (9) 

where $ is an exact solution to the eigenvalue equation (Eq. (3)). We add and subtract 
mass terms to obtain 



E 



M 



E <w+E /W+E v jk (tk-tj)- ( E ] m*c 

j j j<k 

If we define the mass parameter M* according to 



$).(10) 



M*c 2 = 



E a J • + E ^o m o° 2 + E - £j 



^ M 



(11) 



then the associated matrix element for the modified relative contribution vanishes. The 
idea here is that the contribution of the relative problem now appears in the mass 
parameter associated with the center of mass contribution to the eigenvalue. 



2.3. Eigenvalue relation for a single configuration 

With this definition of the mass parameter the eigenvalue simplifies to 



E = 

We can rotate to obtain 
E = /<&' 

where we used 



E 



M 



-CKi • C 



* „2 



f M 



E^ $/ )\/( M * c2 ) 2 + c2 l p l 2 ' 

3 



$' ~ e <P ' R/R . 



(12) 



(13) 



(14) 
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These are the basic results that we reported in [35] • However, our interpretation now is 
different. The mass parameter M* can allow us to connect the center of mass problem 
for the interacting Dirac particles with the center of mass problem of a set of equivalent 
non-interacting Dirac particles that have a different rest mass. 

2.4- Composite masses 

In general the mass parameter M* for the composite will in general depend on P since 
the Dirac Hamiltonian is in general not covariant (it has not yet been shown that M* 
should be invariant for a covariant model, but we think it likely). Also, in a covarient 
model we would expect the energy-momentum relation to be 



E = ^(M c 2 ) + c 2 |P| 2 (15) 

where Mo is the composite rest mass. In light of this, we could correct the non-covariant 
Dirac model by requiring 



E— ft 



$'W(M*c 2 ) 2 + c 2 \P\ 2 J (M c 2 ) + c 2 \P\ 2 (16) 



In this case, we can recover the nonrelativistic limit 



IP 



2 



E ^ M ^ + kk + - (17) 

3. Finite Basis State Model 

We would like to expand our description to include a finite set of states, and in the 

process we would like a formulation in which the state mass impacts the kinematics. 
We begin by assuming a finite basis state model of the form 

* = E c i% ( 18 ) 

j 

We make use of the variational method to obtain the matrix equation 

Ec = He, (19) 

where the matrix H has individual matrix elements given by 

H kl = (^ k \H\^). (20) 
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The vector c made up of the different expansion coefficients is 



c = 



c 2 

V ; / 

This coefficient vector c should not be confused with the speed of light c. 



(21) 



3.1. Diagonal matrix elements and mass parameters 

To arrange for dynamics with the mass parameter matched to the state, we focus first 
on the diagonal matrix elements and write 



H kk = (<£> fc |# R |$ fc > + ($ k \H T \Q k ), 
where the center of mass and relative parts are 

(<S> k \Hn\<f> k ) = (*, 



(22) 



(23) 



($ fc |# r |$ fe > = ($ 



E «i-c*i+E v"/- 2 • E ^(^-^)- 1 E ffii } m 

j<k 



3 3 3<k \ 3 

If we require the relative contribution to vanish (as discussed above), then the state- 
dependent mass parameter \ /,;'; is consistent with 



* 2 



Mlc 



E a 3 ■ cjT 3 + E ftrnj-c 2 + E - £j) 



Off-diagonal matrix elements 
The off-diagonal matrix elements may be written as 



(25) 



ki 



E a J • ^ 



+ E/3. 



,. m . c 2 + ^ - ^.) 



(26) 



where the first term is associated with the center of mass, and the second term is 
associated with the relative problem. 
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The finite basis eigenvalue relations for a composite in free space can be written in the 
form 



Eck 



$1 



£^*0V« c2 ) 2 + c2 l pl2 



i^k 



It will be convenient to write the off-diagonal matrix element in this case as 

H ki = oi kl ■ (cP) + V kh 
where V k i is the coupling matrix element from the relative part of the problem 

E <*j ■ c *j + E h m i c2 + E v 5k{i k - £,) ®i 

j<k 



Via = (*k r 

3 3 

and where the vectors oL k \ are defined by 



OLkl = ($ 



x -, rrij 
? ~M° Cj 

3 



(27) 



(28) 



(29) 



(30) 



We consider in the Appendix the nonrelativistic reduction of a representative ct^i • cP 
matrix element. 



3.4- Rest frame eigenvalue equations 

The eigenvalue equations above are unfamiliar and moderately complicated. It is useful 
to consider them in the rest frame; in this case, we obtain 



Ec k 



$k )M* k c 2 



Ck + E VkiQ. 

lyik 



(31) 



For such a finite basis state model we would expect the basis state masses to appear; 
consequently we may write 



M k c 2 = ($ 



This allows us to write 



^ M 

3 



)M* k c 2 . 



(32) 



E 



C-2 



M x c 2 V 12 V 13 
V 21 M 2 c 2 V 23 
V 31 V 32 M 3 c 2 



V 



/ 



V ; J 



(33) 



It will be convenient to think of this as the basic unperturbed problem 
Ec = H c. 



(34) 
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We can expand the square root terms to lowest order to obtain an eigenvalue equation 
relevant to the low-momentum case; we write 



Ec k 



M k c 2 + 



2Mi 



+ 



i^k 



a k i ■ (cP) + V h 



ci, 



(35) 



where we have approximated the Dirac model as covariant in the estimates of the 
diagonal matrix elements. In matrix notation, this might be written as 

12 

Ec = Mc a + M -1 ^- + a- (cP) • c, (36) 

where 



M 



and where 



Mc 2 


-t-M" 


-l|P| 

2 


( Mi 











M 2 











M 3 


V ; 









"12 


"13 


"21 





"23 


«31 


«32 






(37) 



V : 



(38) 



: / 



4. A Hamiltonian for Nuclei in a Lattice 

There has been discussion over the years as to develop a suitable formalism that would 
be capable of systematically addressing the anomalies of interest in condensed matter 
nuclear science. It was proposed in [37] that one begin with a fundamental Hamiltonian 
based on nucleons and electrons, and then reduce the model for applications by first 
building up nuclei from nucleons, then solving for electronic degrees of freedom in a 
Born-Oppenheimer picture, and finally focusing on the vibrational problem. We imagine 
a similar approach here, only instead of starting with a nonrelativistic fundamental 
Hamiltonian we use a relativistic one. 

The separation of center of mass and relative degrees of freedom is straightforward in 
the nonrelativistic version of the problem, which is why we focused on it in [37]. But we 
see in the discussion above that it is possible to separate the center of mass and relative 
Hamiltonians for a many-particle Dirac model, even in the context of a finite basis 
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approximation. This separation allows us to extend the earlier program systematically 
to a relativistic formulation (including now a new relativistic coupling between the 
nuclear motion and internal nuclear degrees of freedom) based on an underlying many- 
particle Dirac model. 

4-1. Hamiltonian for electrons and nucleons 

We begin with a formal model based on many-particle Dirac Hamiltonians for the 
electrons and nucleons 



In the first term in brackets we find a relativistic nucleon Hamiltonian including 
mass, velocity, and potential terms (including strong force interactions as well as 
electromagnetic interactions). Nuclear models of this kind can be found in the literature 
[3H],[39j. In the second term in brackets we find a relativistic electron Hamiltonian also 
including mass, velocity and potential terms (in this case electromagnetic interactions). 
Electronic models of this kind provide the foundation for relativistic electron band 
calculations; one can find a discussion of this model in [ID]. Although there is no reason 
to believe that a relativistic description for the electrons is required for the problems 
of interest, somehow it seems appropriate to maintain the same level of description in 
the fundamental Hamiltonian. Finally, the last term includes electron-nucleon potential 
terms, which are electromagnetic here. A further augmentation of the model to include 
weak interaction physics is possible, but we will not pursue it here. 

4-2. Reduction of the nucleon Hamiltonian 

The developments presented in the previous sections allow for a systematic reduction 
for the nucleon Hamiltonian in the form 




+ Y, a 3' C Pi + fa™^ + E Vj£ (r* " Tj ) 



electrons 



+ E^ n (r fc -r,).(39) 



j j<k 



E «i • cPi + 0jM jC 2 + E VJT(r k - rj ) 



- nucleons 





(40) 
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where the sum over nucleons is now replaced with a sum over nuclei, and where the 
different matrices associated with the nuclear finite basis expansion are selected to be 
appropriate for the nucleus indexed by I. 

The excited state energies in the rest frame appear as the eigenvalues of the mass 
term Mc 2 . The lowest-order contribution to the kinetic energy of the nucleus as a 
composite Dirac particle is included in |P| 2 /2M, which allows for the energy of a basis 
state to impact the kinetic energy appropriately. Finally, there is a new term a • (cP) 
that describes a relativistic effect in which the nucleus center of mass momentum (which 
will subsequently be part of the lattice vibrations) is coupled to transitions between the 
different basis states. The summation over / here indicates a sum over the different 
nuclei in the lattice, so that there will be separate mass matrices and lattice coupling 
terms for each nuclei. The condensed matter Hamiltonian that results is 



H 



E 



Mc 2 + ^ + a-(cP) 



J nuclei 



electrons 



E^(r fc -r,-)(41) 



j j<k 

This represents a generalization of the conventional starting place for condensed matter 
physics calculations. 



4-3. Born-Oppenheimer approximation 

With nucleons replaced by nuclei, the resulting model is very similar to the standard 
condensed matter model, and we can similarly make use of the Born-Oppenheimer 
approximation to obtain a potential model for the nuclei 

IPI 2 



H 



E 



Mc + ^T + a ' (cP) 



+ ^V J N k N (R k -R j ). 

j<k 



(42) 



This Hamiltonian is made up of individual mass, kinetic energy, and lattice coupling 
terms for each nucleus individually, and augmented now with effective potential 
interactions (electromagnetic plus electronic) between the nuclei. This model provides 
a generalization of the usual lattice Hamiltonian to include nuclear mass effects and 
lattice coupling with the nuclei. 
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4-4- Reduction to a standard lattice Hamiltonian 

In the event that the effects associated with nuclear excitation are weak, then the new 
terms in the model can be dispensed with; if we we assume 

a, -> 0, (43) 

then we recover a model that is essentially the standard condensed matter lattice 
Hamiltonian: 



i 



IP 



2 



H = E Mc 2 + L_L + £ Vj N k N (R k - R,). (44) 



j<k 

In this approximation there is no direct coupling between the nuclear excited states and 
lattice vibrations. If the nuclei are in ground states, we could replace the mass matrices 
by the mass eigenvalues, which completes the reduction: 



i 



IP 



h = E M ^ 2 + 1 df +E^r(R*-R*)- ( 45 ) 

1 j<k 



2M 

5. Discussion and Conclusions 

This study was motivated by our interest in deriving a lattice Hamiltonian from a 
relativistic starting Hamiltonian in which we could study the effect of the different 
configuration masses on the lattice dynamics, with the goal of developing a systematic 
description of the anomalies associated with condensed matter nuclear science. The 
derivation of such a Hamiltonian from the many-particle Dirac model in particular is 
in general problematic due to difficulties in the separation of relative and center of 
mass degrees of freedom for the relativistic problem. We recently obtained a weaker 
result that seemed to suggest a systematic way of separating center of mass and relative 
contributions to the energy, which prompted us to use it to isolate nuclear center of mass 
and nuclear terms in a lattice Hamiltonian. We obtained terms for the mass energy and 
kinetic energy that had been expected. In addition we obtained an additional coupling 
term due to the ex. ■ cP terms in the Dirac model; this we had not anticipated. The 
presence of this new term is not an artifact of our derivation, as we found subsequently 
that it can also be obtained in a straightforward manner using perturbation theory. 

The origin of this effect is that the nuclear states of the composite nucleus transform 
under a boost in the many-particle Dirac model, which implies a mixing with other 
states. In the case of constant P the eigenvalue problem is seeking to create a version 
of the boosted wavefunction out of the rest frame states. However, with the composite 
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momentum is dynamical (as occurs in a lattice), the model tries to develop boosted 
wavefunctions for a composite momentum that keeps changing magnitude and direction, 
which requires in this picture a dynamical admixture of different rest frame states. 

In the end, we obtain a Hamiltonian that describes lattice dynamics and nuclear 
excitation that is derived consistently from an underlying relativistic Hamiltonian 
(the many-particle Dirac model). This model appears to be very closely related to 
models for two-level systems coupled to an oscillator that we have investigated over the 
years [28], [29], [31)], [31] in connection with the excess heat effect in the Fleischmann-Pons 
experiment. 

Appendix A. Nonrelativistic limit for the Transition Operator 

Transitions in this model are described by the off-diagonal matrix element 



OCfi CP = ($/ 



> — -OLj ■ CP 

y m 3 



*<). (a.i; 



We are interested in developing a nonrelativistic approximation for this operator which 
may be useful for understanding the coupling. 

Appendix A.I. Expansion of the wavefunction 

We assume that the solution to the relative problem can be expanded in the form 

$ = $+++... + $_++... + $+_+... + $++-... + ■••, (A.2) 

where the first term involves large component channels for all nucleons, where the second 
involves a small component channel for the first particle, and so forth. 

Appendix A.2. Channel equations 

For the first term, the relative eigenvalue problem in the rest frame results in 



j<k 



$+++... = cr 1 ■ («Ti)$_ ++ ... + er 2 • (c7T 2 )$+_ H 



+ <r 3 • (c7r 3 )$++-... + - • • (A.3) 
For the second term we have 



E 



Mc 2 + 2mic 2 



E 

j<k 



CTi • (c7Ti)$ +++ ... + <T 2 • (C7T 2 )$- 
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+ cr 3 • (c7r 3 )$_ + _... + • • • 
from which we approximate the channel wavefunction as 



$_++... 



E - Mc 2 + 2m lC 2 -J2V jk (£ k - 



16 
(A.4) 

L 

(<r 1 -c7Ti)$ +++ ...(A.5) 



Approximate transition matrix element 

To proceed we expand the transition matrix element in terms of the different pieces 



4^ M J 



mi 
M 



$/(+ + + •••) 



>i ■ cP) 



$,(- + + •••) 



+ (*/(- + + •••) 



(o-i • cP) 



$,(+ + + •••) 



+ 



(A.6) 



We expect the large component to dominate, so we keep terms with a single small 
component and approximate according to 



x - mj -5 
> -rrOCj-cP 



+/$/(+++••• 



> r cP) 



£-Mc 2 +2m lC 2 -]TV- fc (£ fc -^) 
£-Mc 2 +2m lC 2 -^ V$ fc (£ fc -^ 



(o-l-OTi) 



$*(+++•••) 



+ 



$*(+++•••) 

(A.7) 



Taylor series expansion 

A Taylor series expansion yields 



£-Mc 2 +2m lC 2 -]>>, fc (£ fc -^) 



1 



2m ic 2 



1 - 



2m 1 c 2 
E - Mc 2 1 



i E-Mc 2 1 

i+ 2m lC 2 2mT^^ €j 



2mic 2 2mic 2 



(A. 
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Leading-order contribution 

To evaluate the leading-order term in the expansion of the transition matrix element, 
we require the identity 

(a-cr)(b-cr) = a • b + icr ■ a x b. (A. 9) 

This allows us to write 



(<T • CP)(<7 • C7r) + (<T • C7t)((T ■ cP) = 2c 2 7T • P + i(T ■ ( P X 7T + 7T X P 
= 2C 2 7T-P 

Consequently, the leading-order term vanishes 



(A.10) 



( cr ^ cP )^rT2( 'i- c ^) 
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$<(+ + + •••) 



(ATI) 



since 



(A.12) 



Fzrsi next-order term 



There are two next-order terms; the first of these is 



r — ^ m 7 / _ / i . - . _E — Mc 2 . „ . 



(£ - Mc 2 
2Mc 2 



<M+ + + • 



E 



$;(+ + + •••) 
Tf.-.P 



$,(+ + +•••) 



7T ■ 

Note that this term contains the sum J2j — which we may rewrite as 



3 m. 



= E*d^- 



(A.13) 



(A.14) 



since X)j ^"j — 0, where m a „ might be an appropriate average of the different masses. 
Since the proton mass is 938.27 MeV/c 2 and the neutron mass is 939.56 MeV/c 2 , 
the individual nucleon masses are little different from the average nucleon mass. 
Consequently, we expect this term to be small. 
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The other next-order term is 



3 



>3 ■ ^) 



+(*/(+ + + •••) 



tr, • cP) 



2Mc 2 



$/(+ + + •••) Efo-cP) 



1 



$,(+ + + •••) 
(oycP) $,(+++•••) 



.(A.15) 
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